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We describe an entanglement purification protocol to gen- 
erate maximally entangled states with high efficiencies from 
two-mode squeezed states or from mixed Gaussian continu- 
ous entangled states. The protocol relies on a local quantum 
non-demolition measurement of the total excitation number 
of several continuous variable entangled pairs. We propose an 
optical scheme to do this kind of measurement using cavity 
enhanced cross-Kerr interactions. 

PACS numbers: 03.67.Hk, 42.50.-p, 03.65. Bz 
Quantum communication, such as quantum key distri- 
bution and quantum teleportation, is hampered by the 
difficulty to generate maximally entangled states between 
distant nodes [1]. Due to loss and decoherence, in reality 
we can only generate partially entangled states between 
distant sides [2]. Entanglement purification techniques 
are needed to concentrate maximally entangled states 
from partially entangled states [3,4]. For qubit systems, 
efficient entanglement purification protocols have been 
found [3-5] . But none of these purification schemes have 
been realized experimentally due to the great difficulty 
of performing repeated collective operations in realistic 
quantum communication systems. Thus, it is of interest 
to consider purification of continuous variable entangle- 
ment. The nonlocal Gaussian continuous variable en- 
tangled states (i.e., states whose Wigner functions are 
Gaussians) can be easily generated by transmitting two- 
mode squeezed light, and this kind of entanglement has 
been demonstrated in the recent experiment of continu- 
ous variable teleportation [6] . As the first choice for per- 
forming continuous entanglement purification, one would 
consider direct extensions of the purification schemes for 
qubit systems. But till now in these extensions, no entan- 
glement increase has been found for Gaussian continuous 
entangled states [7]. Thus, the discussion should be ex- 
tended to a larger class of operations to purify continuous 
entangled states. Braunstein et al. [|| have proposed a 
simple error correction scheme for continuous variables. 
However, it is not clear whether it can be used for pu- 
rification. In [9] a protocol to increase the entanglement 
for the special case of pure two-mode squeezed states 
has been proposed, which is based on conditional photon 
number subtraction; the efficiency, however, seems to be 
an obstacle for its practical realization. 



In this paper, we present an entanglement purification 
scheme with the following properties: (i) For pure states 
it reaches the maximal allowed efficiency in the asymp- 
totic limit (when the number of pairs of modes goes to in- 
finity); (ii) It can be readily extended to distill maximally 
entangled states from a relevant class of mixed Gaussian 
states which result from losses in the light transmission. 
Furthermore, we propose and analyze a scheme to im- 
plement this protocol experimentally using high finesse 
cavities and cross-Kerr nonlinearities. Our purification 
protocol generates maximally entangled states in finite 
dimensional Hilbert spaces. The entanglement in the 
continuous partially entangled state is transformed to the 
maximally entangled state with a high efficiency. We be- 
gin the paper by describing the entanglement purification 
protocol for pure two mode squeezed states, and then ex- 
tend the protocol to include mixed Gaussian continuous 
states, and last describe the physical implementation of 
the purification protocol. 

First assume that we have generated m entangled pairs 
Ai,Bi (i = 1, 2, • • • m) between two distant sides A and 
B. Each pair of modes Ai,Bi are prepared in the two 
mode squeezed state 1^)^.^., which in the number basis 
has the form 
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where A = tanh (r) , and r is the squeezing parameter 
For and only for a pure state, the entanglement is 
uniquely quantified by the von Neumann entropy of the 
reduced density operator of its one-component. The en- 
tanglement of the state (1) is thus given by E (\^)ab) = 



cosh (r) log (cosh (r)) 



sink 2 



r)log(sinh 2 (r)). The 



joint state l^)^.^-} °f the m entangled pairs is simply 
the product of all the \^) A , which can be rewritten as 
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where (AiBi} is abbreviation of the symbol Ai,B\, 
A2,B2, ■■■ and A m ,B m , and the normalized state 
\i) {AiBi} IS defined as 
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The function & l) in Eq. (2) and (3) is given by 



Am) _ \3-t-m i)i rp Q C oncentrate entanglement of these 
m entangled pairs, we perform a QND measurement of 
the total excitation number n Al + tla 2 + • • • + nA m on 
the A side (we will describe later how to implement 
this measurement experimentally). The QND measure- 
ment projects the state \^)^a b } onto a two-party max- 
imally entangled state \j)^ Ai Bi} wrtn probability pj = 
(l — A 2 )™ A 2j /j m ' ) .The entanglement of the outcome 

state \j) (At B t } is g iven b y E (\j)(A t B t }) = lo g(/j m) )- 
The quantity Tj = E (\j) (AiBi} ) /E (|* ) A . B J defines 
the entanglement increase ratio, and if Tj > I, we get a 
more entangled state. Even with a small number m, the 
probability getting a more entangled state is quite high. 
It can be easily proven that if m goes to infinity, with unit 
probability we would get a maximally entangled state 
with entanglement mE (\^)ab) • This ensures that this 
method is optimal in this limit, analogous to the purifica- 
tion protocol presented in [3] for the qubit case. For any 
finite number of entangled pairs, the present purification 
protocol is more efficient than that in [3], since it takes 
advantage of the special relations between the coefficients 
in the two-mode squeezed state. 

An interesting feature of this entanglement purification 
protocol is that for any measurement outcome j ^ 0, we 
always get a useful maximally entangled state in some 
finite Hilbert space, though the entanglement of the out- 
come state Ij)^.^.} does not necessarily exceed that of 
the original state l^)^^. if j is small. It is also inter- 
esting to note that a small alternation of this scheme 
provides a useful method for preparing GHZ-like states 
in high dimensional Hilbert spaces The key point 

is that the modes Bi need not be at the same side in 
the protocol. Assume we have two entangled pairs B. A\ 
and A 2 , C distributed at three sides B, A, C, with each 
pair being prepared in the state (1). Then a local QND 
measurement of the modes A\ , A 2 at the A side with the 
outcome j ^ generates a three-party GHZ state in the 
(j + 1 )-dimensional Hilbert space. Obviously, if we have 
m entangled pairs, we can generate a (m + I )-party GHZ 
state using this method. 

In reality, the light transmission will be unavoidably 
subjected to loss, and then we will not start from an 
ideal two mode squeezed state, but instead from a mixed 
state described by the following master equation 



(j+m-l)! 
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where p is the density operator of the m entangled pairs 
with p(0) = \^)( Ai Bi} the ideal two mode squeezed 
state, and the effective Hamiltonian 



H cS = -i 
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In Eqs. (4) and (5), a ai denotes the annihilation operator 
of the mode cti (a — A or £?), and we have assumed that 
the damping rates r\ A and r\B are the same for all the m 
entangled pairs based on symmetry considerations, but 
r\A and T)b may be different to each other. 

In many practical cases, it is reasonable to assume 
that the light transmission noise is small. Let r denote 
the transmission time, then tja t and t\bt are small fac- 
tors. In the language of quantum trajectories |f(J, to the 
first order of tjat and t]bt, the final state of the m en- 
tangled pairs is either |* (0) ) (A . B } oc e^ iH " aT \^} {Az b z } 
with no quantum jumps occurred, or |\]A a *)^ ca b } ^ 
y/rjaTa ai \*&)(a-b-} w ith a jump occurred in the on chan- 
nel (a = A,B and i = 1, 2, • • • , m). The final density op- 
erator is a mixture of all these possible states. To purify 
entanglement from the mixed state, we perform QND 
measurements of the total excitation number on both 
sides A and B, and the measurement results are denoted 
by j A and js , respectively. We then compare ja and js 

through classical communication, and keep the outcome 

(i) (i) 

state if and only if ja — 3b- Let P A and P B ' denote 
the projections onto the eigenspaces of the corresponding 

m ra 

total number operators ^ a A a Ai and a B ,aSi with 

i=l * * i=l 

eigenvalue j, respectively. It is easy to show that 
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the outcome state is the maximally 
entangled state |j) ( AiBi } w ith entanglement log^/j" 1 '^. 
The probability to get the state \j)( AiBi \ is now given by 
'1 - A 2 )™ \ 2 3 f( m) e -^A+VB)T]\ It should be noted 
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So if j A = j B = j, 
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that the projection operators P A ' P B ' cannot eliminate 
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the states obtained from the initial state l^) 



(AiBi} 



quantum jump on each side A and B. The total prob- 
ability for occur of this kind of quantum jumps is pro- 
portional to m 2 n 2 riAT]BT 2 . So the condition for small 
transmission noise requires m 2 ri 2 r]Ai]BT 2 <C I, where 
n = sinh 2 (r) is the mean photon for a single mode. 

In the purification for mixed entanglement, we need 
classical communication (CC) to confirm that the mea- 
surement outcomes of the two sides are the same, and 
during this CC, we implicitly assume that the storage 
noise for the modes is negligible. In fact, that the stor- 
age noise is much smaller than the transmission noise is 
a common assumption taken in all the entanglement pu- 
rification schemes which need the help of repeated CCs 
[4,5]. If we also make this assumption for continuous 
variable systems, there exists another simple configura- 
tion for the purification protocol to work. We put the 
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generation setup for two-mode squeezed states on the 
A side. After state generation, we keep the modes Ai 
on side A with a very small storage loss rate tja, and 
at the same time the modes Bi are transmitted to the 
distant side B with a loss rate t\b ^S> t]a- We call this 
a configuration with an asymmetric transmission noise. 
In this configuration, the purification protocol is exactly 
the same as that described in the above paragraph. We 
note that the component in the final mixed density op- 
erator which is kept by the projection PY'P^' should 
be subjected to the same times of quantum jumps on 
each side A and B. We want this component to be a 
maximally entangled state. This requires that the to- 
tal probability for sides A and B to subject to the same 
nonzero times of quantum jumps should be very small. 
This total probability is always smaller than utjat, de- 
spite how large the damping rate r/s is. So the working 
condition of the purification protocol in the asymmetric 
transmission noise configuration is given by nrjAT <C 1. 
The loss rate r]B can be large. The probability to get 
the maximally entangled state \j)r^ iBi \ is still given by 
V 3 = (1 - A 2 )™ A^He-^A+^r/ * 

For continuous variable systems the assumption of 
storage with a very small loss rate is typically unreal- 
istic. If this is the case, then we can use the following 
simple method to circumvent the storage problem. Note 
that the purpose to distill maximally entangled states 
is to directly apply them in some quantum communica- 
tion protocols, such as in quantum cryptography or in 
quantum teleportation. So we can modify the above pu- 
rification protocol by the following procedure: right after 
the state generation, we take a QND measurement of the 
total excitation number on side A and get a measurement 
result ja- Then we do not store the outcome state on side 
A, but immediately use it (e.g., perform the correspond- 
ing measurement as required by a quantum cryptography 
protocol jl^] ) . During this process, the modes Bi are be- 
ing sent to the distant side B, and when they arrive, we 
take another QND measurement of the total excitation 
number of the modes Bi and get a outcome js- The 
resulting state on side B can be directly used (for quan- 
tum cryptography for instance) if 3 a = 3b , and discarded 
otherwise. By this method, we formally get maximally 
entangled states through posterior confirmation, and at 
the same time we need not store the modes on both sides. 

To experimentally implement the above purification 
scheme, we need first generate Gaussian continuous en- 
tangled states between two distant sides, and then per- 
form a local QND measurement of the total excitation 
number of several entangled pairs. Here we propose a 
promising experimental scheme, which uses high finesse 
optical cavity to carry continuous entangled states and 
cavity enhanced cross Kerr interactions to realize the lo- 
cal QND measurement. It is possible to generate Gaus- 
sian continuous entangled states between two distant cav- 



ities [[13| . We can transmit and then couple the two out- 
put lights of the nondegenerate optical parametric ampli- 
fier to distant high finesse cavities. The steady state of 
the cavities is just a Gaussian continuous entangled state 
described by the solution of Eq. (4) after taking into ac- 
count of the propagation loss [14]. The difficult part is to 
perform a QND measurement of the total photon number 
contained in several local cavities. We use the setup de- 
picted in Fig. 1 to attain this goal. (For convenience, we 
use the two-cavity measurement as an example to illus- 
trate the method. Extension of the measurement method 
to multi-cavity cases is straightforward.) 
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FIG. 1. Schematic experimental setup to measure the total 
photon number m + 112 contained in the cavities I and II. The 
cavities I and II, each with a small damping rate k and with 
a cross Kerr medium inside, are put respectively in a bigger 
ring cavity. The ring cavities with the damping rate 7 are 
used to enhance the cross Kerr interactions. A strong cotin- 
uous coherent driving light bn (t) is incident on the first ring 
cavity, whose output b i is directed to the second ring cavity. 
The output b 2 (t) of the second ring cavity is continuously 
observed through a homodyne detection. 

The measurement model depicted in Fig. 1 is an ex- 
ample of the cascaded quantum system 1 10 . The inci- 



dent light bn can be expressed as bn = b n + .9^/7, where 
9\fl (3 is a large dimensionless factor) is a constant driv- 
ing field, and b ix is the standard vacuum white noise, 
satisfying (b£ (t)b' a (V)) = and (b' a (f) $ (7)) = 

S (t — t ^ . The Hamiltonian for the Kerr medium is as- 
sumed to be Hi = %x n ib\bi, (i = 1 or 2) , where 6, is the 
annihilation operator for the ring cavity mode, and x is 
the cross-phase modulation coefficient. The self-phase 
modulation can be made much smaller than the cross 
phase modulation with some resonance conditions for the 
Kerr medium, and thus is negligible [15,16]. In the frame 
rotating at the optical frequencies, the Langevin equa- 
tions describing the dynamics in the two ring cavities 
have the form 
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with the boundary conditions (see Fig. 1) bi 2 = b Q \ = 
hi + .9\/7 + \/7^i and b o2 = b i2 + y/jb 2 . In the realistic 
case 7 3> x( n i) ; (* = 1, 2), we can adiabatically elimi- 
nate the cavity modes bi, and express the final output 
b Q 2 of the second ring cavity as an operator function of 
the observable n\ + n 2 . The experimentally measured 
quantity is the integration of the homodyne photon cur- 
rent over the measurement time T. Choosing the phase 
of the driving field so that g = i\g\, the measured ob- 
servable corresponds to the operator 



the driving field, imbalance between the two ring cavi- 
ties, light absorption of the Kerr media and the mirrors, 
self phase modulation effects, light transmission loss be- 
tween the ring cavities, and inefficiency of the detectors. 
To realize a QND measurement, the imperfections should 
be small enough. We have deduced quantitative require- 
ments for all the imperfections listed above [18]. With 
the parameters given in the above paragraph, all these 
requirements can be met experimentally. 
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where 



and &t ) 
rT , 



satisfying 

1, is defined by br — ^7= f$ b a (t) dt. Equa- 
tion (8) assumes 7 ^> \ (m) and e~ lT <§; 1. There are 
two different contributions in Eq. (8). The first term rep- 
resents the signal, which is proportional to n\ 4- n 2l and 
the second term is the vacuum noise. The distinguisha- 

If 

we effec- 

1.1 </ - \- 

tively perform a measurement of n\ 



bility of this measurement is given by Sn — , , . 
Sn < 1, i.e., if the measuring time T > r ..,~U ; 

- n 2 ; and if T is also 



smaller than 



the photon loss in the cavities I and 



II during the measurement is negligible. So the setup 
gives an effective QND measurement of the total photon 
number operator n\ + n 2 under the condition 



7 



< T < 



1 



(9) 



This condition seems to be feasible with the present tech- 
nology. For example, if we assume the cross Kerr in- 
teraction is provided by the resonantly enhanced Kerr 
nonlinearity as considered and demonstrated in [15,16], 
the Kerr coefficient x/^ n ~ 0.1MHz would be obtain- 
able [17]. We can choose the decay rates k/2it ~ AM Hz 
and 7/27T ~ 100MHz, and let the dimensionless factor 
g ~ 100 (for a cavity with cross area S ~ 0.5 x 10 _4 cm 2 , 
g ~ 100 corresponds a coherent driving light with in- 
tensity about 40mWcm~ 2 ). The mean photon number 
(rii) = (712) = sinh 2 (r) ~ 1.4 for a practical squeezing 
parameter r ~ 1.0. With the above parameters, Eq. (9) 
can be easily satisfied if we choose the measuring time 
T ~ 8ns. More favorable values for the parameters are 
certainly possible. 

To bring the above proposal into a real experiment, 
there are several imperfect effects which should be con- 
sidered. These imperfections include phase instability of 
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